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We explore the cosmological constraints on the parameter Wdm of the dark matter 
barotropic equation of state (EoS) to investigate the "warmness" of the dark matter fluid. 
The model is composed by the dark matter and dark energy fluids in addition to the radia- 
tion and baryon components. We constrain the values of w<$m using the latest cosmological 
observations that measure the expansion history of the Universe. When u>d m is estimated 
together with the parameter Wd c of the barotropic EoS of dark energy we found that the 
cosmological data favor a value of u>dm = 0.006 ± 0.001, suggesting a warm dark matter, 
and Wdc = —1-11 ± 0.03 that corresponds to a phantom dark energy, instead of favoring a 
cold dark matter and a cosmological constant (u>d m = 0, Wd c = — 1)- When Wd m is estimated 
alone but assuming w dc = -1, -1.1, -0.9, we found w dm = 0.009 ± 0.002, 0.006 ± 0.002, 
0.012 ± 0.002 respectively, where the errors are at 3er (99.73%), i.e., Wdm > with at least 
99.73% of confidence level. When (u>dm, ^dmo) are constrained together, the best fit to data 
corresponds to (w,j m = 0.005 ± 0.001, QdmO = 0.223 ± 0.008) and with the assumption of 
Wdc = —1-1 instead of a cosmological constant (i.e., u>dc = — !)• With these results we found 
evidence of if dm > suggesting a warm dark matter, independent of the assumed value for 
Wdc, but where values Wd c < — 1 are preferred by the observations instead of the cosmolog- 
ical constant. These constraints on Wd m are consistent with perturbative analyses done in 
previous works. 

PACS numbers: 04.20.-q, 04.70.Bw, 04.90.+C 

Keywords: Warm dark matter, cosmological observations, constraints 



* avelino@fisica.ugto.mx 
t norman.cruz@usach.cl 

* ulises@ifm.umich.mx 



2 



I. INTRODUCTION 

The astrophysical evidence for the existence of Dark Matter (DM) is well based on observations 
from the scales of galaxies, clusters and the universe itself, in the framework of the standard 
cosmological model. 

Despite the fact that the cosmological scenario where cosmological parameters fit a dark matter 
mainly non-baryonic and cold, a great debate has currently opened about the possibility that 
Warm Dark Matter becomes a better candidate to understand the recent investigations. For a wide 
discussion based in the new results in the area see pj|. A summary of astrophysical constraints on 
dark matter is present in [2] 

Let us summarize some of the difficulties of the Cold Dark Matter (CDM) model. At galactic 
scales, N-body simulations of cosmological structures with CDM have predicted that the dark 
matter halos surrounding galaxies must present radial profiles of the mass density and velocity 
dispersion with a central cusp in which the value of the logarithmic slope is under discussion (see 

In the case of the missing satellites problem [4|], [5] exists a discrepancy in the cold dark matter 
model between the predicted numbers of satellite galaxies inside the galactic halo for the Milky 
Way and lower number observed. This problem has been undertaken assuming a warm dark matter 
component in various works j^j. 

A recent investigation of N-body simulations in a warm dark matter scenario, which took into 
account the new dwarf spheroidal galaxies discovered in the Sloan Digital Sky Survey (SDSS) Q], 
derived lower limits on the dark matter particle mass [8]. 

In recent investigations, the measuring of the dark matter equation of state (EoS) has been 
carried out using different approaches. Following a suggestion given in [9], where the method 
combines kinematic and gravitational lensing data, the dark matter EoS was measured in 1(| 
using galaxy clusters which present gravitational lensing effects. The result of this work indicates 
that the measured EoS for dark matter is consistent with the standard pressureless cold dark 
matter at la level. Nevertheless, lensing analysis in clusters such as Coma and CL0024 shows a 
trend to prefer an exotic EoS for the dark matter, i.e., w ~ —1/3. 

Models of the dark matter component described by a fluid with non-zero effective pressure has 
been studied in some astrophysical scenarios. At galactic level, an EoS with anisotropic pressures 
has been investigated in [111 ] in order to explain flat rotation curves. A polytropic dark matter 
halo fits very well a number of elliptical galaxies, improving or at least giving similar results to the 
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velocity dispersion profile compared to a stars-only model 
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Explorations of the EoS for dark matter at cosmological level have been carried out in various 
frameworks. In [13], a constant EoS for dark matter is studied from the study of the power spec- 
trum, assuming a cosmological constant as the dark energy fluid and a flat universe. The bounds 
obtained for Wd m were —1.50 x 10 -6 < Wd m < 1-13 x 10 -6 if there is no entropy production and 
—8.78 x 10~ 3 < Wd m < 1-86 x 10~ 3 if the adiabatic sound speed vanishes. Phenomenologically, 
EoS for both dark fluids have been studied in []jj] using WMAP+BAO+i^o observations by syn- 
chronizing the model with the ACDM model at the present time. The dark matter component 
behaves like radiation at very early times and at the present time Wdm = 0.0005. 

In the case of unified dark matter models, where a single matter component is assumed to source 
the acceleration and structure formation [15], the initial phase is described by a cold dark matter 
so that the fitting with cosmological data leads to a late phase with negative Wdm very close to a 
cosmological constant or phantom matter. 

Our aim in this work is to study the EoS of the dark matter component allowing a non zero 
value for Wd m from the beginning and then to undertake a constraining of its value using the latest 
observations that measure the expansion history of the universe. In what follows, we will assume 
a barotropic EoS for this component. Of course, this assumption is rather restrictive because in 
approaches based in the nature of particles constituting the dark matter fluid is expected to have 
a Wdm varying with the cosmological time. Such is the case, for example, for dark matter Bose- 
Einstein condensation [if]]. Nevertheless, if the dark matter fluid is modeled, in the non-relativistic 
approximation, as a non-degenerated ideal Maxwell-Boltzmann gas, a barotropic EoS is obtained 
with Wd m = constant [itI ] . 

The present paper is organized as follows. In Section II, we briefly outline the basic equations 
of evolution of the model. In Section III, the parameters of the model are constrained using cosmo- 
logical data from type la supernovae, CMBR, baryon acoustic oscillations, the Hubble expansion 
rate and the age of the universe. Finally, in Section IV, we discuss and conclude our results. 



II. THE COSMOLOGICAL MODEL 

We study a cosmological model composed by four fluids: radiation, baryons, dark matter and 
dark energy. We assume a barotropic equation of state (EoS) for dark matter (dm) and energy 
(de) fluids, pi = u>, • pi, with i = dm, de, respectively, pi corresponds to the density of the fluid 
and pi to its pressure. We are interested in studying the cosmological prediction for the EoS of the 
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dark matter, in particular, for the magnitude of w,j m . 

We assume a spatially flat Friedmann-Robertson- Walker (FRW) cosmology. The Friedmann 
constraint and conservation equations for the radiation, baryonic, dark matter and dark energy 
fluids are given respectively as 



H =— — (Pr + Pb + Pdc + Pdm) (1) 

=p\ + AHp T (2) 

=p\ + 3H Ph (3) 

=Pd m + 3Hp dm (l + w dm ) (4) 

=p dc + 3Hp dc (l + w de ), (5) 



where H is the Hubble parameter and the dot over p\ stands for the derivative with respect to the 
cosmic time. The conservation equations ©-(Jl]) have the respective solutions in terms of the scale 
factor a 



I \ PrO I \ PrO / \ PdmO / s PdeO / c x 

Pl (a) = -^, Pb(a) = ^ ; Pdm(a) = a3(1+wdm) , Pde ( a ) = (6) 

where the subscript zero at pio indicates the present-day values of the respective matter-energy 
densities. Inserting the expression on the Friedmann constraint ([T]), and dividing by the Hubble 



constant Hq, it becomes 



p 2 ( \ = H2 ( a ) _ SttG /pro PbQ PdmO , PdcO 

{a) ~ Hi ~ 3H$ V a 4 a 3 a^+^nO + a 3(i+™ dc ) 



(7) 



We define the dimensionless parameter densities as = Pm I Pent > where p® rit is the critical den- 
sity evaluated today defined as p® rit = 3Hq /(8ttG). With this definition, the Friedmann equation 
(|7|) obtains the form 



-p2f \ _ ^rO ^bO ^dmO ^deO ,„\ 

W ~ a 4 a 3 a 3 ( 1+w 'd m ) + a 3(i+^ dc ) ' W 
or, using the relation between the scale factor and the redshift u z" given by a = 1/(1 + z), we 
rewrite the dimensionless eq. ([8]) in terms of the redshift as 



E\z) = Q r0 (l + zf + Q h0 (l + zf + fWoU + zf 1+w ^ + Q de0 (l + z ) 3 ( 1 +^ e ). (9) 
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Setting E{z = 0) = 1 we have the constraint equation, 

^deO = 1 — (^r0 + ^b0 + ^dmo)- (10) 

III. COSMOLOGICAL CONSTRAINTS 

To constrain the value of w dm using cosmological data, to compute their confidence intervals and 
to calculate their best estimated values, we use the following cosmological observations described 
below measuring the expansion history of the Universe. 

To perform the numerical calculations, it was used for the baryonic and radiation (photons and 
relativistic neutrinos) components the values of f^o = 0.0458 [3] and f2 r o = 0.0000766 respectively, 
where the later value is computed from the expression {^(J 



Q l0 = n 7 „(i + 0.2271iV s 



cfi 



(11) 



where N e g = 3.04 is the number of standard neutrino species 



19, 



2J and fi 7 o = 2.469 x 10~ 5 /i 



51,-2 



corresponds to the present-day photon density parameter for a temperature of T cm b = 2.725 K 
where h is the dimensionless Hubble constant h = #o/(100 km/s-Mpc). 



1. Type la Supernovae 



We use the type la supernovae (SNe la) of the "Union2.1" data set (2012) from the Supernova 



Cosmology Project (SCP) composed of 580 SNe la 
flat FRW Universe is defined as 



22]. The luminosity distance di in a spatially 



dL(z,w dr 



c(l + z) 



dz' 



(12) 



H Jo E(z',w dm ) 

where "c" corresponds to the speed of light in units of km/sec. The theoretical distance moduli ^ 
for the k-th supernova at a distance Zk is given by 



^\z,w dra ) = 5 log 



dh{z, w dm ) 
Mpc 



+ 25 



(13) 



So, the x 2 function for the SNe la test is defined as 
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Best estimates for (wdm,Wd e ) 



Data set 


Wdxn 




Xmin 


Xd.o.i. 


SNe la 




n nnfi +0133 


-1.003i°;^ 


562.23 


0.97 


(H,l A ,z m ) CMB 




0.004 ±0.001 


+0.057 
-0.053 


1.11 


1.11 


H(z) 




07+ 069 

u - uu ' -0.058 


-1 197+ 014 

J" 15 " -0.13 


8.05 


0.73 


SNe + CMB + BAO + 


H(z) 


0.006 ±0.001 


-1.115 ±0.033 


578.84 


0.97 



TABLE I. Best estimated values for (mdm,Wdc)- See figures Q] and [2] for the confidence intervals. We find 
that the cosmological data used in the present work favor a non-vanishing magnitude, positive value for 
Wd m suggesting a warm dark matter, in addition to wj c < — 1 indicating a phantom dark energy. In order 
to compare these results with the ACDM model we computed the value of the % 2 function evaluated at 
(wMm = 0,u>d c = —1) using the same four cosmological data sets (SNe + CMB + BAO + H(z)) together, 
finding a value of Xacdm = 740.5, that is clearly greater than Xmin = 578.8 obtained in the present work 
for Wdm — 0.006, Wdc = —1.115, indicating that the ACDM model fits not too well the cosmological data 
compared with the latter values. It was assumed f2t>o = 0.0458, O r o = 0.0000758, OdmO = 0.23 and Hq = 73.8 
km/(sMpc). The errors correspond to 68.3% of confidence level (la). 

Best estimates for Wdm 





2 




Wdm 


w dc Xmin 


Xd.o.f. 


0.009 ±0.002 


-1 591.57 


0.99 


0.006 ±0.002 


-1.1 579.02 


0.97 


0.012 ±0.002 


-0.9 628.21 


1.05 



TABLE II. Best estimated values for Wdm when it is assumed the different values of Wdc = —1,-1.1, —0.9 
for the dark energy. In the three cases it is found a non- vanishing positive value for Wdm- We find also that 
the best fit to data is for the case when it is assumed Wdc = —1-1, i-e., it has the smallest value of Xmin 
compared with the other two cases. It was used the joint SNe + CMB + BAO + H(z) data sets. The errors 
are at 99.73% of confidence level (3c). See figure |3] for the likelihood functions. 



XSNel^dm, H ) = 2^1 I (14) 

k=l \ °« / 

where is the observed distance moduli of the k-th supernova, with a standard deviation of 
in its measurement, and n = 580. It was used a constant prior distribution function for Hq to 
marginalize it (i.e., it is not assumed any particular value of Hq) because Hq is a nuisance parameter 
in the SNe la test. 
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Confidence Intervals, SNe Union2.1 



Confidence Intervals, H(z) 




0.000 0.002 0.004 0.006 0.008 




FIG. 1. Confidence intervals for (w,jm, Wdc)- The upper left panel corresponds to the use of the SNe la data 
set release (2012) "Union 2.1" of the SCP 22j, through the minimization of the x 2 function (fl4|) . The upper 
right panel corresponds to the use of Hubble parameter data at different redshifts using the x 2 function 
(|3"U|) . The lower left panel corresponds to the use of three observational data (TZ, Ia,z*) given by WMAP-7y 
19], through the x 2 function (|22|) . And the lower right panel corresponds to the use of the total % 2 function 
([31]) that contains the four type of cosmological observations together SNe + CMB + BAO + H(z). The 
best estimated values for (w^m, w^c) of each panel are indicated with the red point and the magnitudes are 
shown in table HI It is assumed a spatially flat FRW universe and for the baryon, dark matter, radiation 
parameter densities and the Hubble constant it was assumed the values of Obo = 0.0458, fidmo = 0.23, 
n iQ = 0.0000766 



n 



19, 21] and H = 73.8 km/s-Mpc [l8| respectively. The contour plots correspond to 68.3% 



{la), 95.4% (2a) and 99.73% (3a) of confidence level. 
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FIG. 2. Confidence intervals (CI) all together for (wdm, Wdc) calculated with the different cosmological data 
sets (cf. figure [I] for CI separately for each cosmological data set). The CI labeled by "R-CMB" and "do. 275 
BAO" come from the use of the shift parameter 1Z and the distance ratio d z at z = 0.275 of BAO computed 
through the x 2 functions defined at (fl7|) and (|29|) respectively. The right panel corresponds to a zoom in 
of the left panel in order to show the CI coming from the use of the (7£, Ia>z*) CMB distance priors and 
the joint SNe+CMB+BAO+i/(z). The best estimated values are shown in table HI The interval regions 
corresponds to 68.3% (la), 95.4% (2a) and 99.73% (3a) of confidence level. 
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FIG. 3. Likelihood functions for u>d m when it is assumed the values of Wdc = —1, —1-1, —0.9 for the EoS of 
dark energy. See table |TT] for the best estimated values of Wdm- It is assumed a spatially flat FRW universe 
and for the baryon, dark matter, radiation parameter densities and the Hubble constant it was assumed the 
values of il h0 = 0.0458, Sl dm0 = 0.23, O r0 = 0.0000766 and H Q = 73.8 km/s-Mpc respectively. 
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Best estimates for (f«dm,^dmo) 



Data set 




2 2 
Xmin A-d.o.f . 


SNe la 


0.004 ±0.27 0.229±g:J| -1 
-0.103 0.316 -1.1 
0.177 0.138 -0.9 


562.22 0.972 
562.20 0.972 
562.275 0.972 


[K,Ia,z*) cmb 


-0.0009 ±0.003 0.183±g:gjg -1 
0.001 0.206 -1.1 
-0.007 0.153 -0.9 


2.558 2.558 
0.04 0.04 
10.90 10.90 


H(z) 


n2 -|t-+0.226 nn4 +0.076 i 
U.ZlO_ 0212 U.114_ 050 1 

0.085 0.176 -1.1 
0.426 0.057 -0.9 


8.271 0.751 
8.12 0.738 
8.55 0.777 


SNe + CMB + BAO + H(z) 


0.005 ±0.002 0.204 ±0.008 -1 
0.005 ± 0.001 0.223 ± 0.008 -1.1 
0.005 ± 0.002 0.184 ± 0.007 -0.9 


582.11 0.978 
578.27 0.971 
601.41 1.010 



TABLE III. Best estimated values of the parameter density of dark matter fidmo and the Wdm of the EoS 
of dark matter (pdm = ^dro ' Pdm)- The hrst column shows the cosmological data sets used to compute the 
best estimates shown in second and third columns. The fourth column indicates the assumed value for u>d c - 
The fifth and sixth columns correspond to the minimum value of the x 2 function, Xmin' ana - X 2 by degrees 
of freedom, Xd f respectively. The latter is defined as \\ f = Xmin/( n ~ P)> where n is the number of data 
and p the number of free parameters (in this case p — 2). The computed values come from the minimization 
of the x 2 functions defined in {HJ, J22]), (00]) and J31} respectively. The fourth row (SNe + CMB + BAO 
+ H(z)) encloses the information of all the cosmological observations used in the present work to constrain 
the values of (f2dm0> Wdm)- Notice that Wdm has a positive value, favoring a warm instead of a cold dark 
matter. See figures l4l to [TOl for the confidence intervals. 

2. Cosmic Microwave Background Radiation 



We use the WMAP 7- years distance priors shown in table 9 of [19(, composed of the shift 
parameter 1Z, the acoustic scale I a and the redshift of decoupling z*. 
The shift parameter 1Z is defined as 



n 



-(1 + z*)Da(z*) 



(15) 



where Q m o corresponds to the total present pressureless matter in the Universe, i.e. f2 m o = Qbo + 
^dmO) and Da is the proper angular diameter distance given by 

dz' 



D A (z) 



{l + z)H J E(z',w d 



(16) 
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Confidence Intervals SNe Union2.1 Confidence Intervals, (R, l A ,z,) CMB 




0.2 0.4 0.6 0.8 1.0 0-16 0-17 018 0-19 0.20 0.21 0.22 

fl-dmO ^dmO 



Confidence Intervals, H(z) SNe + CMB + BAO + H(z) 




0.0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.18 0.19 0.20 0.21 0.22 0.23 



FIG. 4. Confidence intervals for the present-day value of the parameter density of dark matter fidmO 
versus the Wd m of the EoS of dark matter, and assuming w,j c = —1 for the dark energy. The upper 
left, right, lower left and right panels correspond to the use of SNe, (11,1a, z#) of CMB, H(z) and the 
joint SNe+CMB+BAO+i? (z) data sets respectively. The best estimated values are indicated with the red 
point and the magnitudes are shown in table Mil It is assumed a spatially flat FRW universe and for the 
baryon, dark matter, radiation parameter densities and the Hubble constant it was assumed the values of 
ftbo = 0.0458, ft dm0 = 0.23, Q r0 = 0.0000766 and H = 73.8 km/s-Mpc respectively. The interval regions 
correspond to 1, 2 and 3cr of confidence level. 



for a spatially flat Universe. With TZ we can denned a \ 2 function as 

XR-CMB(^dm,#o) = f — a ^ 0bS ) ( 17 ) 

where 7£ bs = 1.725 is the "observed" value of the shift parameter and on = 0.018 the standard 
deviation of the measurement (cf . table 9 of [3] ) . 
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The acoustic scale I a is defined as 



Ia = (1 + z*) 



ttDa(z*) 
r s {z*) 



(18) 



where r s (z*) corresponds to the comoving sound horizon at the decoupling epoch of photons, z*, 
given by 



r s (z) 



da 



(19) 



VSJo a 2 F(a)^l + (30bo/4^ 7 o)a 

where as mentioned above, we use f2 7 o = 2.469 x 10~ 5 /i -2 as the present-day photon energy density 
parameter, and fi&o = 0.02255/i~ 2 as the baryonic matter component, as reported by Komatsu et 



al. 2011 

Sugiyama 



where 



19(] . We compute the theoretical value of z* from the fitting formula proposed by Hu and 



23] 



= 1048 [l + 0.00124(fi bo /i 



2\-0.7381 



[i + 9l (n m0 h 2 y*] 



9i 



0.0783(fi b o^ ) 



2n,-0.238 



0.560 



l + 39.5(O b o/i 2 ) - 763 ' 92 l + 21.1(fi b0 / i 2 ) 1 -8 1 ' 
The x 2 function using the three distance priors (Ia,TZ, -z*) is defined as 

3 



(20) 



(21) 



XCMEiKim, #o) = (%i ~ di)(C 



ij(xj dj) 



(22) 



where Xi = (Ia,7Z, z*) are the theoretical values predicted by the model and di = (I a = 302.09, 1Z = 
1.725, z* = 1091.3) are the observed ones. For Hq it was assumed the latest reported value of 
Hq = 73.8 km/s-Mpc [18]. The is the inverse covariance matrix with entries |]J 



C 



-i 



2.305 29.698 -1.333 
29.698 6825.27 -113.180 
1.333 113.180 3.414 



(23) 



3. Baryon Acoustic Oscillations 

We use the baryon acoustic oscillation (BAO) data from the SDSS 7- years release |2 !]. expressed 
in terms of the distance ratio d z at z = 0.275 defined as 

r s {zd) 



d, 



0.275 



D v (0.275) 



(24) 
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Cosmological Constraints Cosmological Constraints (Zoom 1) 




FIG. 5. Confidence intervals (CI) for (fidmOj Wdm), calculated with the different cosmological data sets (see 
figure 2]). It is assumed a value of Wdc = —1 for the parameter of EoS of dark energy. The right panel 
corresponds to a "zoom in" of the left panel in order to show the tiny CI that come from the use of the 
(TZ, Ia, z*) CMB distance priors through the Xcmb function and labeled as "CMB". The CI from the 
total x 2 function (|31 [) are even smaller than those of the CMB and so they are shown in figure [HI The best 
estimated values are shown in table Hill The interval regions corresponds to 68.3% (la), 95.4% (2a) and 
99.73% (3<t) of confidence level. 



where z^ is the redshift at the baryon drag epoch computed from the fitting formula 25] 

(^rno/i 2 ) - 251 



z d = 1291- 



l + h^oh 2 )^ 



1 + 0.659(ft m0 /i 2 ) - 828 
h = 0.313(^ m0 /i 2 )" - 419 [1 + 0.607(^0 /i 2 ) ' 674 ] 



b 2 = 0.238(Q m0 /r) 



2\0.223 



(25) 

(26) 
(27) 



For a flat Universe, Dy(z) is defined as 



D v (z) 



dz' 

W7) 



Hz) 



1/3 



(28) 



It contains the information of the visual distortion of a spherical object due the non-Euclidianity 
of the FRW spacetime. 

The value do b 2 s 75 contains the information of the other two pivots, drj.2 and usually used 

by other authors, with a precision of 0.04% 
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Cosmological Constraints (Zoom 2) 




I Ll I u I ^3 u I I ^ J I I 0^ I I 

0.16 0.18 0.20 0.22 0.24 



FIG. 6. Confidence intervals for (^dmO: ^dm)- This figure corresponds to a "zoom in" of the figures [4] and 
[5] to show the CI that come from the use of all the observational data sets together through the total x 2 
function defined in ([31]) and labeled as "SNe + CMB + BAO + H(z)" . The best estimated values computed 
with the total \ 2 function are u>dm = 0.005, rWiO = 0.204 (see table [TTT| ) . The intervals regions correspond 
to 68.3% (la), 95.4% (2a) and 99.73% (3<r) of confidence level. Notice that w dm > with 95% confidence 
level, suggesting a warm dark matter. 



The x function for BAO is defined as 



2 , „ , , , _ /' d 0.275 - do b 275 X ^ 



xWwdm, H ) = ^ r^_^ j (29 ) 
where dn^ s 7 rr = 0.139 is the observed value and ad = 0.0037 the standard deviation of the measure- 



ment [24]. For Hq it was assumed the latest reported value of Ho = 73.8 km/s-Mpc 18]. 



4- Hubble expansion rate 

For the Hubble parameter, we use the 13 available data, 11 data come from the table 2 of Stern 
et al. (2010) [26] and the two following data come from Gaztanaga et al. 2010 Q: H(z = 0.24) = 
79.69 ± 2.32 and H(z = 0.43) = 86.45 ± 3.27 km/s/Mpc. For the present value of the Hubble 
parameter, we take the value reported by Riess et al 2011 [18]: H(z = 0) = Hq = 73.8 ± 2.4 
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FIG. 7. Confidence intervals for (^dmOj^dm) when it is assumed the value of w^c = —1.1 for the parameter 
of EoS of dark energy, i.e., a phantom dark energy. See table Hill for the values of the best estimates. The 
interval regions corresponds to 68.3% (la), 95.4% (2a) and 99.73% (3c) of confidence level. 



km/s/Mpc. The \ 2 function is defined as 

X K (w dm , H ) = 2^ \ — I (30) 

where H(zi) is the theoretical value predicted by the model and Hf hs is the observed value with 
its standard deviation au- 

Finally, with the \ 2 functions defined above we construct the total \ 2 function given by 

X 2 = XsNe + XcMB+xlAO+XH- (31) 
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FIG. 8. Confidence intervals for (^dmOj^dm) when it is assumed the value of w& c = —1.1 for the parameter 
of EoS of dark energy, i.e., a phantom dark energy. See also figure [71 Table Hill shows the values of the best 
estimates for this case. The central and right panels correspond to a zoom in of the left panel. The interval 
regions corresponds to 68.3% (la), 95.4% (2a) and 99.73% (3c) of confidence level. 



Age of the Universe 



Age (Gyr) 


Wdm 




^dmO 


Table Point 


13.23 ±0.06 


0.006 


-1.115 


0.23 
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A 


12.97 ±0.012 


0.009 


-1 


0.23 
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B 


13.19 ±0.015 


0.007 


-1.1 


0.23 


El 


C 


12.7 ±0.01 


0.013 


-0.9 


0.23 


El 


D 


13.37 ±0.14 


0.005 


-1 


0.204 


ED 


E 


13.31 ±0.12 


0.005 


-1.1 


0.223 


urn 


F 


13.4 ±0.12 


0.005 


-0.9 


0.184 


ED 


G 



TABLE IV. Age of the Universe given in gigayears (first column) when it is assumed certain values for 
{wdm, Wde, ^dmo) shown in the 2nd to 4th columns and that comes from the best estimates of Wdm shown 
in tables H HlIII f fifth column). The last column indicates the letters used in figure [11] to label those points. 
Hq is assumed to be 73.8 km/s-Mpc. 



We minimize this function with respect to the set of parameters (u>dm, ^de), (Vdm>^dmo) an d 
Wdm alone, to compute their best estimated values and confidence intervals or likelihood functions. 
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FIG. 9. Confidence intervals for (fidmOj w d m ) when it is assumed the value of w^e = —0.9 for the parameter 
of EoS of dark energy. See table Hill for the values of the best estimates. The interval regions corresponds 
to 68.3% (la), 95.4% (2er) and 99.73% (3a) of confidence level. 



5. The age of the Universe 

Using the fact that H = a /a, we can rewrite the eq. (JSJ) as an ordinary differential equation 
(ODE) for the scale factor a in terms of the cosmic time as 



d a qtt /^rO . ^bO . ^deO ^dmO n / on x 

— - pH a\ —t- H r + —— r + — j— r = 0, (32) 

dt V a 4 a 3 a 3(i+w dc ) a 3(l+w dm ) y ' 

where (3 = 1.022729 x 1CT 3 is introduced to give the units of time in gigayears (Gyr) when the 
value of the Hubble constant is given in units of km/(s-Mpc). For the conversion of units, we use 



the values of 1 year = 31558149.8 seconds (a sidereal year)[29j] and 1 Mpc = 3.0856776 x 10 19 km 
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FIG. 10. Confidence intervals for (fidmO> WMm) when it is assumed the value of w,j c = —0.9 for the parameter 
of EoS of dark energy. See also figure [9l Table IIIII shows the values of the best estimates for this case. The 
central and right panels correspond to a zoom in of the left one. The interval regions corresponds to 68.3% 
(la-), 95.4% (2a) and 99.73% (3a) of confidence level. 



29|], so $ = (31558149.8 x 10 9 )/3.0856776. 

We solve numerically the ODE (|32|) with the initial condition a(t = 0) = 



30J and compute 



the value ttoday of the age of the Universe through the condition a(t to da y ) = 1- 

Evaluating the numerical solution of the ODE (|32p at the best estimates and assuming the 



values of H = 73.8 ± 2.4 



], Q r0 = 0.0000758, n h0 = 0.0458 ± 0.0016 Q we find an age of the 
Universe. See table [IV] and figure [TTJ From the oldest globular clusters the age of the Universe is 
constrained to 12.9 ± 2.9 Gyr [28j. 



IV. DISCUSSION AND CONCLUSIONS 



We explored the constraints on the value of the parameter Wdm of the barotropic EoS of the 
dark matter to investigate the "warmness" of the dark matter fluid. The model is composed 
by the dark matter and dark energy fluids in addition to the radiation and baryon components. 
We constrained the value of w^ m using the SNe la "Union 2.1" of the SCP data set, the three 
observational (TZ, Ia, z*) data from the CMB given by WMAP-7y, the distance ratio d z at z = 0.275 
of BAO and the Hubble parameter data at different redshifts. 

We calculated the best estimated values for the pair of parameters (wdm, Wd e ), (wdm, ^dmo) and 
also Wdm alone, where Wd e and ^dmo are the parameter of the barotropic EoS of dark energy and 
the present-day value of the density parameter of dark matter respectively. 

When Wdm is estimated together with Wd e we found that the cosmological data prefer the 
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FIG. 11. Age of the Universe given in gigayears (Gyr) as a function of Wdm- The black long dashed, the solid 
red and the blue short dashed lines correspond to assume the values of Wdc = — 1, — 1-1, —0.9 respectively. 
The right panel corresponds to a zoom in of the left one, where the points locate the inferred value of 
the age of the Universe when the eq. (1521) is evaluated at the best estimated values for Wdm (see section 
IIII 51 and table [TV)) . The letters that label the points correspond to the values of (wd mi Age, u>d e > ^dmo) 
where: A = (0.006,13.23,-1.115,0.23), B = (0.009,12.97,-1,0.23), C = (0.007,13.19,-1.1,0.23), 
D = (0.013,12.7,-0.9,0.23), E = (0.005,13.37,-1,0.204), F = (0.005,13.31,-1.1,0.223) and G = 
(0.005, 13.4, —0.9, 0.184). The shaded area corresponds to the consistent region for the age of the Universe 
estimated from the oldest globular clusters (Age= 12.9 ± 2.9 Gyr 



m) 



value of Wdm = 0.006 ± 0.001, suggesting a warm dark matter, and Wde = — 1.11 ± 0.03 that 
corresponds to a phantom dark energy, instead a cold dark matter and a cosmological constant 
(wdm = 0,Wde = — !)• See table U and figures [TJ and [2j 

In order to study the dependence of the estimated value for Wdm with respect to the value of Wde 
of the dark energy, we computed the best estimate of Wdm as the only free parameter but assuming 
three different values of u>d e = —1,-1.1,-0.9. We found the values of u>d m = 0.009 ± 0.002, 
0.006 ± 0.002, 0.012 ± 0.002 when it is assumed the values of Wd e = —1,-1.1,-0.9 respectively, 
where the errors were computed at 3er (99.73%), so, we found that w^m > with at least 99.73% 
of confidence level. Additionally, from these three cases, the assumption of Wd e = —1.1 is the case 
that allows to fit better the model to data compared with the other two cases (see table HI] and 
figure [3|) . 

When Wd m is constrained together with ^dmo we found that the best fit to data is for (u'dm = 
0.005 ± 0.001, SldmO = 0.223 ± 0.008) and with the assumption of u>de = —1-1, instead of a 
cosmological constant (i.e., Wd e = — 1)- We found also interesting to notice that the best estimated 
value of Wdm using all the combined data sets give the same value of Wdm = 0.005 independent of 
the assumed value for Wd e , where the three cases were Wd e = — 1, — 1-1, —0.9 (see the three rows at 
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the bottom of table Ullj) . 

In all cases the best fit to data, measured through the Xd o f magnitude, of the cosmological 
observations separately or all together (the joint SNe + CMB + BAO + H(z) data) correspond to 
the case when it is assumed w& e = — 1.1 (phantom dark energy) instead of a cosmological constant 
(w Ae = -1) or w dc = -0.9. 

For the age of the Universe, we found a consistent value for the age when it is evaluated at the 
best estimated values for w^m, except for the case when it is assumed w dc = —0.9. See table ITvl 
and figure [TTJ 

On the other hand, Muller [3] and more recently Calabrese et al. 31] investigated the con- 



straints on Wdm at perturbative level comparing with the large scale structure data and CMB 
anisotropies. They found the constraints —0.008 < Wdm < 0.0018 and —0.0133 < Wd m < 0.0082 
respectively. We find that our results are comparable and consistent with these ones. 

In summary, we found an evidence of a non- vanishing value if dm- From the cosmological obser- 
vations we found constraints on the values of u>d m around 0.005 < w^m < 0.01 suggesting a warm 
dark matter, independent of assumed value for w Ae , but where a value Wd e < — 1 is preferred by the 
observations instead of the ACDM model. Our constraints on w,im are consistent with perturbative 
analysis done in previous works. 
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